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Abstract 

The interface tension in the three-dimensional Ising model in the 
low temperature phase is investigated by means of the Monte Carlo 
method. Together with other physically relevant quantities it is ob- 
tained from a calculation of time-slice correlation functions in a cylin- 
drical geometry. The results at three different values of the temper- 
ature are compared with the predictions from a semiclassical approx- 
imation in the framework of renormalized (/)^ theory in three dimen- 
sions, and are in good agreement with them. 
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1 Introduction 



In various systems in statistical mechanics or in nature at temperatures T 
below a certain critical temperature interfaces (domain walls) may be 
present, separating coexisting phases. The interface tension r is the free 
energy per unit area of interfaces. It is an interesting quantity for various 
reasons. For numerical simulations of field theories in a phase with broken 
discrete symmetry finite volume effects are often dominated by tunneling 
[Q]. The size of these finite volume effects is intimately related to the inter- 
face tension. Secondly, the interface tension is relevant for first order phase 
transitions, where it determines the nucleation rate. 

For many binary systems the interface tension r has been investigated 
experimentally as a function of the temperature T ||, |, || . As T increases 
towards Tc the reduced interface tension 

where k is Boltzmann's constant, vanishes according to the scaling law 

a ~ aot^ , (2) 

where 

T - T 

(3) 



t 



T 

c 

and (To is the critical amplitude of the interface tension. Widom's scaling law 

/i = 2i/, (4) 

relates the universal critical exponent fi to the critical exponent of the cor- 
relation length ^: 

^r^^ot-'. (5) 

The experimental results H, ^, Q for the critical index /i, 

/i = 1.26 ±0.01, (6) 

and for the critical exponent z/ p|, ||, |l^ , as well as the prediction from the 
renormalization group [|r^, 

= 0.630 ± 0.002 , (7) 



2 



are consistent with Widom's law. 

Other quantities, which are universal according to the scaling hypothesis 
[ P^ , [T^, are the dimensionless products of critical amplitudes 



where and C,o are the critical correlation length amplitudes in the high- 
temperature and low-temperature phases, respectively. 

On the theoretical side we dispose of two models in the relevant uni- 
versality class, the Ising model and the 0^ field theory in three Euclidean 
dimensions. Early field theoretic calculations [|14| , |15| of -R_ (and R^) in 
the framework of the e-expansion, where d = 4 — e is the number of di- 
mensions of space, were in disagreement with the experimental results. The 
discrepancy could, however, be resolved by a calculation in renormalized field 
theory directly in three dimensions [1^. This calculation was performed by 



means of a semiclassical saddle-point expansion to one-loop order around the 
kink-solution. The semiclassical method yields an expression for the energy 
splitting between the two "vacua" due to tunneling in a finite volume, which 
involves the interface tension. For some quantities related to the energy split- 
ting and the interface tension it yields definite predictions in terms of the 
renormalized parameters of the model (see below). 

A similar development took place in the case of Monte Carlo simulations 
of the Ising model. First results for the interface tension |]T^ , which were at 



variance with the (real) experiments, have been improved by recent numerical 
studies which made the discrepancy vanish. 



A systematic numerical study of the interface tension and related objects, 
which would allow a comparison with the semiclassical approximation, was. 



however, completed only recently for the four-dimensional Ising model 

In this article we present the results of a Monte Carlo investigation of the 
three-dimensional Ising model. The interface tension is determined from the 
volume dependence of the energy splitting, a method that was introduced 
ST[] . Furthermore, masses and renormalized parameters are calculated, 



m 



which are necessary for a comparison with the semiclassical formulae. 

Let us add a remark about a controversy, which stems from work of Borgs 
and Imbrie [Q. They derive a formula for the asymptotic large volume be- 



haviour of the energy splitting, which disagrees with the semiclassical result. 
Their result, however, is valid only at very low temperatures, where interfaces 
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are rigid. In three dimensions this low temperature region is separated from 
the scahng region in the low temperature phase by the roughening transition 
Above the roughening transition there are massless interface modes and 



22 and the 



the results of ||22| do not apply to the rough phase. 

In four dimensions, however, there is a conflict between 
semiclassics, because of the absence of a roughening transition. The Monte 
Carlo results of are in favour of the semiclassical formula. Nevertheless 



it might be possible that Borgs and Imbrie are right for very large volumes. 
Namely, if one adopts the common assumption that the interface dynamics 
is effectively described by a SOS-model or equivalently a discrete Gaussian 
model, the results of Mack and Gopfert suggest that in the scaling region 
there is a very small mass in the interface dynamics, which only shows up in 
the behaviour at very large volumes. 



2 The three-dimensional Ising model 

We consider the Ising model on a three-dimensional cubic lattice of cylindrical 
geometry. This means that the lattice has a quadratic cross-section of area 
in two directions and extends over T lattice units in the third direction, 
such that L <^T . (T is not to be confused with the physical temperature, 
which does not appear any more from here on.) The Boltzmann factor is 
denoted 

e-^^ , (9) 

with ^ 

^ = -EE</'-0.+A , /3>o, (10) 

X fl=l 

where ft denotes the unit vector in the positive //-direction. The variables (p^ 
are associated with the lattice points and assume values (px = ±1- They obey 
periodic boundary conditions. The transfer matrix is denoted , where H 
is the Hamiltonian in the language of quantum theory. The lowest eigenvalue 
of H is normalized to zero. 

The three-dimensional Ising model in an infinite volume is known to pos- 
sess a critical point at 



/3c = 0.221 654 (6) . (11) 
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For (3 > f3c the Z2-symmetry — —0 is spontaneously broken and the field 
acquires a non-zero vacuum expectation value, e.g. 

(0,) = ±v, v> 0. (12) 

In this phase the spectrum of the transfer matrix is doubly degenerate. In 
particular there are two ground states or "vacua" |0±) with 

(O+|0(x)|O+) =t;, (O-|0(x)|O_) = -V. (13) 

On the other hand it is well known in statistical mechanics that sponta- 
neous symmetry breaking does not occur in a finite volume. Let us consider 
a situation in which L is finite. As a consequence of the Frobenius-Perron 



theorem applied to the transfer matrix there is a unique ground state 
\0s) symmetric under the refiection —(f), and the vacuum expectation 
value of the field vanishes. This means that the degeneracy of the infinite 
volume ground states |0±) is lifted. Separated from the ground state |0s) by 
a small energy splitting Eoa there is an antisymmetric state |0a), and if one 
decomposes these states as 

|0,) = (|0+) + |0_)) |0„) = i= (|0+) - |0_)) , (14) 

then |0+) and |0_) are states which go over into the degenerate vacua in the 
infinite volume limit. 

The energy splitting Eoa is due to tunneling between |0+) and |0_) in a 
finite volume. Its volume dependence was studied in ^]. Their analysis 
is based on a picture of domains which extend over the spatial volume and 
cover certain intervals in time. Neighbouring domains with a different sign of 
the field are separated by domain walls, which can be considered as tunneling 
events. From this picture a prediction about the energy splitting of the form 

Eo„~exp{-aL2} (15) 

is obtained, where a is the interface tension associated with the domain walls. 
One sees that tunneling effects vanish very rapidly with increasing volume. 

The next higher states are the symmetric and antisymmetric one-particle 
states with energies Eu and Eia, respectively. In the limit L oo they 
become degenerate too, and their energy defines the physical mass m. 
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In a numerical simulation the low-lying spectrum of zero-momentum 
states can be obtained from correlations of time-slices 



Let 



Tr e 



-TH 



(16) 



(17) 



be the partition function. Then the vacuum expectation value of the product 
of timeslice field averages is given by 

-tEoa g-(T-t)Eoaj _^ 



= V 



+clo 



-01 



g-tEla _^ ^-{T-t)Ela 
{T~t)Eoa-tEi, " 



^-tEia-(T-t)E^s ^-{T-t)Eia-tEu 



+ 



(18) 



where the matrix elements are defined as 

V={0s\St\0a) 
Coi = {Qs\St\la) Cio = (l^l^'tlOa) 



(19) 

cn = (l,|5t|l„). (20) 

The field expectation value v and the energy splitting Eoa can be extracted 
from the t-dependence of (SoSt), if T and t are chosen suitably large. We 
define the connected time-slice correlation function in a finite volume by 



Here one can use the approximation 

Z =1 + 6"^^°" 

for large enough T. The susceptibihty is given by 

X2 = L'J2{SoSt)c, 



(21) 



(22) 



(23) 
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and the second moment of the correlation function is defined by 



/i,^3L2E(^o^*)cSin^f (sin^^)"'. (24) 

The correlation function of squared time-slices, (SqS^), has an expansion 
similar to (0), from which the energies Eis and i?2s can be obtained. 

Near the critical point, where one is deep enough in the scaling region, 
the critical behaviour can be described in terms of the three-dimensional 0^ 
theory, which is in the same universality class. The quantities, which are 
usually handled in field theory, are defined differently from those, which are 
common in statistical physics. For the definition of the renormalized mass 
mn, the wavefunction renormalization factor Zji and the renormalized field 
expectation value v^i we refer to . The relation to the quantities introduced 
above is in three dimensions given by 



2 6x2 ^ 6/?X2 



VR = ^P/Z~RV. (25) 



Ai2 yU2 
A suitably defined dimensionless renormalized coupling is 

«« = 3^. (26) 

The energy splitting E^a can be calculated in a semiclassical approxima- 
tion, which amounts to a saddle-point expansion around the classical kink 
solution of 0^ theory. The kink interpolates between the two field values 
at the minima of the potential and represents an interface separating re- 
gions with different local mean values of the field. The calculation including 
quadratic fluctuations around the kink solution has been done in |^ for d=4 



dimensions and in [|I6[ for the three-dimensional case. In three dimensions 
the result is 

Eoa = Cexp{-a(L)L^}, (27) 
where a is the interface tension, and 

C = 4£^^,/A,„. (28) 
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The interface tension has a neghgible exponentially small L-dependence and 
the value at L = oo is 

Since m is the inverse of the correlation length in the low temperature phase, 
we obtain an expression for the universal amplitude product 

where mJj is the value of ur at the critical point. From a combination of 
various analytical results in the literature it can be estimated to be 

= 15.1 ±1.3. (31) 

which yields 

/2_ = 0.1024 ±0.0088 (32) 

and 

C/m = 0.49 ± 0.02 (33) 

for (3 ^ (3,. 

3 Monte Carlo calculation 

We performed Monte Carlo simulations of the three-dimensional Ising model 
on lattices with an extension of T = 120 in the third direction. For the 
simulation the algorithm of Swendsen and Wang [^] was employed. Wolff's 



one-cluster version [Q, which usually performs better, was not used, since 
we wanted to take advantage of the improved estimators for four-point func- 
tions. Three different values of the inverse temperature f3 were taken, namely 
/3 = 0.2275, 0.2327 and 0.2391. The size L was varied in steps of 2 in the 
ranges 10-18, 8-14 and 4-10, respectively. From the correlation functions 
of time-slices and squared time-slices the field expectation value v and the 
energies E^a, Eia, Eu and E2s were obtained by means of fits to the t- 
dependence. Furthermore the susceptibility X2 and the second moment /i2 
have been calculated. The results for the energies are displayed in table 1. 
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The table also contains estimates for the physical mass m, which are obtained 
by weighted averages of Eis and Eia on the largest available lattices. For 
the /3-value closest to the critical point the L-dependence of the low-lying 
spectrum is shown in figure 1. 

For each value of P the energy splittings have been fitted according to 

Eoa = Cexp{-aL^}, (34) 

in order to get the interface tension and the prefactor C. The fits are dis- 
played in figure 2 and the numerical results are collected in table 2. 
The values for a are then fitted to the scaling law in both forms 



and 

The results are 
and 



a = ao\l-/3//3,r (35) 

a = ao\l-f3J(3r- (36) 
fx = 1.22(1) , ao = 1.29(4) (37) 



/i = 1.28(1) , ao = 1.64(5) , (38) 

respectively. The fit (^) is drawn in figure 3. The corresponding fits for the 
physical mass, which is equal to the inverse correlation length ^, yield 

u = 0.61(5) , mo = 3.8(6) (39) 

and 

u = 0.64(6) , mo = 4.3(7) . (40) 

The differences are due to the fact that one is not close enough to the critical 
point, but the critical exponents are well consistent with the known values 
from the literature. Forming the universal amplitude ratio we obtain 

R. = 0.09(3) , (41) 

in excellent agreement with the estimate (^). 

From the correlation function of time slices the field theoretic quantities 
have been calculated too. The results are contained in table 3. The second 
moment /i2 has a relatively large error. This is due to the fact that it is 
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sensitive to the correlation function at larger distances, where the signal-to- 
noise ratio is smaller. As a consequence also the estimates for rriR and the 
coupling Ufi are afflicted by large uncertainties. The numerical values for 
rrifi are, however, well consistent with those for m, as expected from series 



expansion results [Q. Furthermore the values for ur agree with the estimate 
( |3TD within errors. 

In order to compare with the semiclassical formulae we formed the ra- 
tios a/m? and C/m from the measured values. The results are shown in 
table 4. Although in the temperature range under consideration a and C 
are varying by factors of 4 and 2, respectively, the dimensionless ratios are 
nearly constant, supporting the semiclassical prediction. The corresponding 
couplings ur, calculated according to and are also shown. As the 
table reveals the couplings obtained in this way are in good agreement with 
each other for each value of (3. They are also consistent with the measured 
couplings as well as with the estimate 



4 Conclusion 

The Monte Carlo simulation of the three-dimensional Ising model on a lattice 
with cylindrical geometry allowed to determine the energy splitting E^a in 
a finite volume with an accuracy sufflcient for a calculation of the interface 
tension a. The scaling law for a and for the inverse correlation length m was 
verified and the resulting amplitude product R_ is in very good agreement 
with an estimate based on a semiclassical approximation in the framework 
of cj)^ theory. The dependence of a and the prefactor C on the renormalized 
mass rriR and on the renormalized dimensionless coupling ur predicted by 
the semiclassical approximation to one-loop order is in agreement with the 
results of the numerical simulation. 

Acknowledgement: The calculations have been performed on the "Lan- 
desvektorrechner" of the RWTH Aachen and on the VAX 9000 of the Uni- 
versity of Miinster. 
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Table 1 



The low-lying energies on lattices of size ■ T . The physical mass m in 
the last row is the weighted average of Eis and Eia on the largest available 
lattices. 







Els 


Ela 


E2s 


m 


(3 = 0.2275 


L =10 


0.04334(8) 


0.42(4) 


0.69(6) 


1.4(4) 




L =12 


0.0217(1) 


0.36(4) 


0.55(3) 


1.0(2) 




L =14 


0.00989(7) 


0.35(4) 


0.455(16) 


0.84(14) 




L =16 


0.00400(7) 


0.38(5) 


0.410(13) 


1.0(2) 




L =18 


0.00151(8) 


0.39(3) 


0.420(12) 


1.0(4) 


0.416(16) 


/3 = 0.2327 


L = 8 


0.0348(1) 


0.54(4) 


0.90(15) 


1.4(2) 




L =10 


0.01069(7) 


0.51(2) 


0.67(2) 


1.2(1) 




L =12 


0.00253(10) 


0.54(2) 


0.606(20) 


1.1(2) 




L =14 


0.00047(20) 


0.59(7) 


0.62(1) 


1.5(5) 


0.619(14) 


(3 = 0.2391 


L = 4 


0.174(1) 


1.05(2) 


5(3) 


1.8(3) 




L = 6 


0.0491(2) 


0.74(7) 


1.35(9) 


2.1(5) 




L = 8 


0.00990(3) 


0.70(6) 


0.86(2) 


1.5(2) 




L =10 


0.00128(6) 


0.79(5) 


0.808(16) 


2.0(6) 


0.806(20) 



Table 2 

The interface tension a and the prefactor C obtained from a fit of the L- 
dependence of Eoa according to (^3|). 



p 


a 


InC 


C 


0.2275 


0.0150(1) 


-1.675(19) 


0.187(4) 


0.2327 


0.0328(2) 


-1.260(13) 


0.284(4) 


0.2391 


0.0572(2) 


-0.955(10) 


0.385(4) 
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Table 3 

Renormalized parameters determined from the numerical simulation. 



/5 


m 


rriR 


vb. 


Zb 


Ur 


0.2275 


0.416(16) 


0.44(05) 


0.262(30) 


0.80(19) 


19(7) 


0.2327 


0.619(14) 


0.62(12) 


0.34(06) 


0.71(27) 


16(9) 


0.2391 


0.806(20) 


0.80(32) 


0.40(16) 


0.65(52) 


15(18) 



Table 4 

The dimensionless ratios a/m? and Cjm obtained in the numerical sim- 
ulation. The columns to the right of these ratios contain the coupling ur 
calculated from the ratio according to the semiclassical expressions ( ^9]) and 
(PHD, respectively. 





/ '1 
a jm 


Ur 


Cjm 


Ur 


0.2275 


0.087(7) 


17.1(1.0) 


0.450(20) 


18.0(1.6) 


0.2327 


0.086(4) 


17.2(0.6) 


0.458(12) 


17.4(0.9) 


0.2391 


0.088(4) 


16.9(0.6) 


0.477(13) 


16.1(0.9) 



Figure Captions 



Fig. 1: The low-lying spectrum at /3 = 0.2275 as a function of the lattice 
size L. 

Fig. 2: The logarithm of the energy splitting E^a versus the lattice cross- 
section for three values of the inverse temperature /?. Also shown are the 
fits according to (|3^). 



Fig. 3: The logarithm of the interface tension a versus ln(l — j3c/ j3). Also 



shown is the fit according to the scaling law (p6D 
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